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ABSTRACT 

In interstellar dust grains, internal processes dissipate rotational kinetic energy. The dissipation 
is accompanied by thermal fluctuations, which transfer energy from the vibrational modes to rota- 
tion. Together, these processes are known as internal relaxation. For the past several years, internal 
relaxation has been thought to give rise to thermal flipping, with profound consequences for grain 
alignment theory. I show that thermal flipping is not possible in the limit that the inertia tensor does 
s ! \ not vary with time. 

Subject headings: ISM: dust 

s : 
<■ 

1. INTRODUCTION 

■ 

■ Two processes act in concert to align grains with the interstellar magnetic field: (1) the grain's principal axis of 
greatest moment of inertia d\ aligns with respect to its angular momentum vector J and (2) J aligns with respect 
i — -i to the magnetic field vector B. Purcell (1979) noted that internal mechanisms for dissipating rotational energy drive 
i-C ■ the grain to its lowest energy state for a given J, namely steady rotation with a\ || J. This occurs on a much shorter 
timcscale than that on which external processes align J relative to B. 
q Purcell (1979) identified two internal dissipation mechanisms. Inelastic dissipation results from the periodic me- 
$_i ■ chanical stresses experienced by a grain that does not rotate steadily about a principal axis. The existence of this 
' process is fairly obvious, but calculating the dissipation rate is a challenging problem (see, e.g., Sharma et al. 2005 
d . and references therein). Purcell (1979) introduced a second, subtle effect, which he termed "Barnett dissipation". 
When a grain does not rotate steadily about a principal axis, the angular velocity vector u) varies periodically in a 
coordinate system attached to the grain. If the grain consists of a paramagnetic material, then the microscopic spins 
(with gyromagnetic ratio j g ) attempt to align with the fictitious "Barnett-equivalent" magnetic field Bbe = <*>/7g- 
As the grain magnetization attempts to follow Bbe, rotational kinetic energy is dissipated. This process is analogous 
to a magnetic resonance experiment, where the dissipated energy is provided instead by the applied radiation field. 
, Purcell (1979) provided a heuristic derivation of the Barnett dissipation rate for oblate grains with dynamic sym- 
£f) • metry. "Dynamic symmetry" refers to the case that Ii = I3, where It are the moments of inertia associated with the 
principal axes a,. Thus, for oblate grains with this symmetry, I\ > 1% = I3. (Henceforth, the term "oblate" shall 
always refer to dynamic, rather than geometric, symmetry.) In this case, Bbe consists of a static component /^g)d\ 
' plus a component Bbe, rot with magnitude Wj_/7 S that rotates in the &2 — S3 plane with angular speed w TO t- Solving 
f*^ the Euler equations yields 

w|l = — C0S7 , (1) 
U' u± = — sin 7 , (2) 

a ; h 

and 

Ah - h) 

^rot = -r- T cos 7 , (3) 

-'l-'2 

where 7 is the (constant) angle between J and a\. Assuming J and 1\ are constant, it is convenient to introduce a 
dimensionless measure of the rotational energy E: 

2/ E 

<7=-^ = l + (r 2 -l)sin 2 7 , (4) 

where the final equality is for oblate grains. Note that, for oblate grains, q ranges from 1 to r% = I1/I2. 
Purcell (1979) argued that the dissipation rate is given by 

d f) =-V X "Bl^ Iot u Iot , (5) 

1 / Bar 

where V is the grain volume and x" is the imaginary component of the magnetic susceptibility. It is worth noting that, 
although this expression (as well as a variant in Lazarian & Draine 1999b) is widely used in grain alignment theory, 
it has not yet been rigorously derived or experimentally verified. Purcell adopted the low-frequency susceptibility 

X" ~ Xow ro t^2 , (6) 
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where xo is the static susceptibility and T 2 is the spin-spin relaxation time. 
Combining equations (JXJ) through |6]) yields 

§ = -r^fa-l)(r a -g) (7) 

with 

TBar " 2 X0 yr 2 J2 ■ (8) 

More realistic approximations for x" (e.g., Draine & Lazarian 1999) yield more complicated expressions for dq/dt, but 
retain the linear dependence on (q — 1) and (r 2 — q) near 9 = 1 and q = r2, respectively. 

In the inverse process of Barnett dissipation, a fluctuation spontaneously transfers some energy from the thermal 
reservoir provided by the grain vibrational modes to the grain rotation. Lazarian & Draine (1997) showed that these 
thermal fluctuations can play an important role in grain alignment. They examined the classic alignment model, 
developed by Purcell (1975, 1979), in which a systematic torque r sys , fixed in grain body coordinates, spins the grain 
to "suprathermal" rotation. Thermal rotation, arising solely from collisions with particles from a gas with temperature 
T gas , is characterized by J ~ Jth = y Jgas (fee is Boltzmann's constant). Suprathermally- rotating grains, with 
J Jth, are impervious to disalignment by random collisions with gas atoms. Thus, J can gradually align with B 
via the Davis-Greenstein (1951) mechanism. Purcell (1979) found that the most important systematic torque results 
from the formation (and subsequent ejection) of H 2 molecules at special sites on the grain surface. 

The distribution of molecule-forming surface sites can change rapidly compared with the Davis-Greenstein alignment 
rate (see, e.g., Lazarian 1995). As a result of this resurfacing, r sys • Si may reverse sign, sometimes spinning the grain 
down to thermal rotation. (In inertial coordinates, T sys is always parallel or anti-parallel to J, when averaged over the 
grain rotation.) These episodes, known as crossovers, were first studied by Spitzer & McGlynn (1979), who concluded 
that thorough disalignment of J relative to B occurs after passage through a small number of crossovers. Lazarian & 
Draine (1997) found that the small disalignment of a% from J during periods of suprathermal rotation (due to thermal 
Barnett fluctuations) limits the minimum value of J during a crossover, thereby limiting the disalignment of J with 
B. 

Although Lazarian & Draine (1997) found that thermal fluctuations during periods of suprathermal rotation may 
aid grain alignment, these same authors soon concluded that thermal fluctuations during periods of slow rotation 
may severely suppress alignment. Lazarian & Draine (1999a) introduced the concepts of thermal flipping and thermal 
trapping. For an oblate grain, internal thermal fluctuations cause the angle 7 between J and Si to vary stochastically 
(see eq. 2]). Whenever 7 crosses n/2 (a "thermal flip"), T sys ■ J changes sign. If flips occur rapidly, then the grain 
can only achieve suprathermal rotation if it reaches, by random walk, a sufficiently high J that the flipping timescale 
(which increases with J) exceeds the spin-up timescale. Grains for which suprathermal rotation is thereby suppressed 
are "thermally trapped" . 

Purcell (1979) considered only the contribution of electron paramagnetism to Barnett dissipation. Lazarian & 
Draine (1999b) showed that nuclear paramagnetism can yield much larger dissipation rates for thermally-rotating 
grains. When including "nuclear relaxation" in their analysis, they found that grains with size up to 1/xm are trapped. 
Thus, the Purcell (1979) scenario of Davis-Greenstein alignment of suprathermally rotating grains appears to fail, 
unless the grains contain superparamagnetic inclusions (Jones & Spitzer 1967). 

Radiative torques (Harwit 1970a, b; Dolginov 1972; Draine & Weingartner 1996, 1997; Weingartner & Draine 2003; 
Hoang & Lazarian 2008; Lazarian & Hoang 2007, 2008), which are not fixed in grain body coordinates, have the 
potential to rapidly align J with B. Weingartner & Draine (2003) found that radiative torques can drive grains into 
various alignment states, some characterized by thermal rotation and some by suprathermal rotation. For the former 
states, the grains were thought to undergo rapid flipping. This result was confirmed by additional calculations in 
Hoang & Lazarian (2008), who also noted that thermally rotating, aligned grains may ultimately reach aligned states 
characterized by suprathermal rotation, due to random gas atom impacts. 

Thermal flipping appears to play a critical role in grain alignment theory, precluding the Purcell scenario (i.e., Davis- 
Greenstein alignment with suprathermal rotation suppressing disalignment) and affecting the aligned grain states in 
the radiative torque scenario. Thus, a quantitative estimate of the flipping rate is needed. This can be accomplished 
with the use of the Langevin and/or Fokker-Planck equations. (Gardiner 2004 provides an excellent introduction to 
stochastic methods.) In §2, I will show that thermal flipping as described by Lazarian & Draine (1999a) is, in fact, 
not possible. 

2. THE LANGEVIN EQUATION FOR INTERNAL RELAXATION 

The Langevin equation is a stochastic differential equation describing the time evolution of the grain rotational 
energy: 

dq = A(q) dt + yj D{q) dw , (9) 

where dw is a Gaussian random variable with variance dt. For Barnett dissipation (in the Purcell 1979 approximation), 
the drift coefficient A(q) is given by the right hand side of equation ([7]). Ideally, the diffusion coefficient D{q) would 
also be derived from the model for Barnett relaxation, but no model has been developed with sufficient detail for this to 
be possible. Instead, D(q) can be determined (to within a constant of integration) by demanding that the probability 
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current S(q) vanish at all q when thermal equilibrium obtains. If f(q)dq is the probability that the dimensionless 
energy lies between q and q + dq, then 



2 dq 



(eq. 5.2.8 in Gardiner 2004). 
Weingartner & Draine (2003) defined the quantity 



8=1-- 
7T 



where 



da 
vr/2 



h{h - hq) + h (h - h) cos 2 a 



h(h - h) + h {h - h) cos 2 a 

, q < h/h 
, q > h/h 



1/2 



h(hq-Il) 
hih-h) 



1/2 



(11) 

(12) 



and showed that the density of energy states is constant in s. This holds for grains with arbitrary h, 7 2 , h- For oblate 
grains (J 2 = h), 

1/2 

■s : 1 - I - — - I . (13) 



ra-1 

Thus, for oblate grains, the thermal equilibrium distribution function is 



where 



/te(<?) oc exp(-kq)— ex exp(-fcg)(r 2 - g) 



J 



-1/2 



(14) 



(15) 



(16) 



2hk B T d 

The thermal equilibrium distribution function is more complicated for grains lacking dynamic symmetry, but still 
depends on k. 

The impossibility of thermal flipping can be simply demonstrated by examining the relaxation at q = r 2 in the 
limit that the dust temperature Td — > 0. In this limit, k — > oo. As Td — > 0, fluctuations cease to contribute to the 
probability current S(q), implying that d(fD)/dq — > 0. This limiting behavior must hold for all q (including q = r 2 ) 
and for any physically realizable probability distribution f(q). Normalization of f(q) requires that any divergence at 
q = r 2 be shallower than f(q) oc (r 2 — q)~ x , unless f(q) — S(q — r 2 ). 

Of course, the contribution of drift to the current, A(q)f(q), must also vanish at q — r 2 . Evidently, it is necessary 
that A(q) falls off linearly or faster with (r 2 — q) near q = r 2 . Note that the Barnett dissipation rate of equation ||7J) 
does satisfy this condition. Equation §8§ suggests that rear — » as Td — > 0, since xo cx TJ" 1 . However, this cannot be 
correct, since A(g = r 2 ) would be undefined rather than zero as Td — > 0. Equation ||5J) does not hold for lower than 
the Curie temperature. For such low temperatures, the material is ferromagnetic, suggesting that rear approaches a 
non-zero constant as Td — > 0. 

Focusing now on the fluctuating term at q = r 2 , 

lim d[f(k,q)D{k,g)] _ Q 

(/c- 1 ,r 2 -g)-»(0,0) C?q 

The limit only exists if it takes the same value for all paths along which (fc _1 , r 2 — q) — > (0, 0). Since there exist paths 
for which (r 2 — q) — ► arbitrarily more rapidly than — > 0, d(f D) / dq may not contain any divergences with respect 
to q = r 2 . Thus, for (r 2 — g) <C 1, fD must either (1) be independent of (r 2 — g) or (2) fall off linearly or faster with 
( r 2 — q)- If condition (1) holds for a particular distribution fi(q), then it will not hold for another distribution / 2 (q) 
having a different dependence on (r 2 — q) near g = r 2 . Thus, condition (2) must generally obtain, implying that D 
must fall off as (r 2 — q) 2 or faster near q — r 2 . This implies that Z? and dD/dq both vanish at g = r 2 . If the diffusion 
coefficient is smooth with respect to k (in the sense that dD/dk exists for all k), then these conditions must be satisfied 
for all k. 

Thus, A, D, and dD/dq all vanish at q = r 2 , making this point a "natural boundary" (see §5.2.1e of Gardiner 2004). 
A system can never reach a natural boundary if it begins at a different point (i.e., with a different value of q). However, 
the system must reach q = r 2 and return to lower q (with a different sign for cos 7) in order for a flip to occur (see 
§2.5.2 of Weingartner & Draine 2003). Consequently, thermal flipping is prohibited. This conclusion does not depend 
on the form of A(q), except that A(q) decreases as (r 2 — q) or faster for q near r 2 . It holds for any type of internal 
relaxation and for grains with or without dynamic symmetry, so long as dD/dk exists for all k. 

If thermal flipping is truly prohibited, then this result must obtain regardless of the choice of variable. Although 
the current S is independent of variable, the two terms composing it, representing drift and diffusion, arc not. When 
transforming variables in stochastic differential equations, the ordinary rules of calculus only apply for linear trans- 
formations. Otherwise, Ito's formula must be used (see §4.3.3 of Gardiner 2004). When the Langevin equation ([9]) is 
transformed to variable y(q), the result is 



dy = 



ai \ d v 1 1 r>( \ d2 y 

A{q) Tq + 2 D ^_ 



dt 



(17) 
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Note the additional contribution to the drift coefficient when the Langevin equation is written in the new variable. 
Ito's formula, along with the relation f(q)dq = f(y)dy, yields 



If y(q) oc (r 2 



q) p (with p =^ 0) for q near r 2 , then 

d 2 y/dq 2 
dy/dq 



oc (r 2 - q)~ 



(18) 



(19) 



Since D(q)f(q) oc (r 2 — g) n with n > 1, the second term in equation (| 18() vanishes at q = r 2 . Thus, if the diffusion 
contribution to the current vanishes at the point q = r 2 for variable q, then it does so for arbitrary variable. 

To illustrate the above arguments in a concrete setting, I will now discuss the diffusion coefficient D(q) for an oblate 
grain and the approximate dissipation rate of equation (|7|) . Setting the probability current equal to zero for thermal 
equilibrium yields 



D(q) = 



1 



/te(?) 



D(1)/ te (1) + 2 / A(q')h E (q')dq' 



(20) 



Upon integrating, 

k 2 r Bal D(q) = [3 



2k(q - l)](r 2 - q) + C(k)(r 2 - q) 1 ' 2 exp(fcg) 



k-^ 2 [3 + 2k{r 2 - l)](ra - <?) 1/2 exp[-fc(r 2 - g)] 



V2, 



^k{r 2 -q) 



exp(x 2 )dx 



(21) 



with 



C(fc) = fc 2 exp(-fc)(r 2 - l)- 1/2 i?(g = 1, k)r Bal - 3exp(-fc)(r 2 - 
+fc~ 1/2 [3 + 2fc(r 2 - 1)] exp(-fcr 2 ) / exp(a; 2 )dx 



1) 



1/2 



For (r 2 - g) < 1, 



r D(g, fc) « C(fc) exp(fcg)(r 2 - q) 1 / 2 + -k 2 (r 2 - l)(r 2 - qf 



(r 2 - q) < 1 



(22) 



(23) 



The term containing C(fc) does not fall off sufficiently quickly with (r 2 — g). Thus, C(k) = identically (for all k, if 
dD/dk exists for all k). The remaining term varies as (r 2 — q) 2 , the shallowest permissible dependence. Note that the 
term containing C(k) satisfies condition (1) following equation (|16|) for the thermal equilibrium distribution /te(<z), 
but not the required general condition (2). 

Given the above general argument prohibiting thermal flipping induced by internal relaxation, one may ask how 
Lazarian & Draine (1999a) concluded that thermal flipping is possible. Their analysis was highly approximate and 
did not employ any diffusion coefficient. Nevertheless, their estimate of the flipping rate agreed well with the detailed 
analysis of Roberge & Ford (1999), which made use of the Barnett relaxation diffusion coefficient calculated by Lazarian 
& Roberge (1997). 

Lazarian & Roberge (1997) solved a modified version of equation (|10p for the diffusion coefficient, in which they 
used the angle 7 rather than q as the variable. They considered oblate grains and the approximate Barnett dissipation 
rate in equation ([7]). Using equation ^ to substitute for q in terms of 7 in equation ([7]), they adopted 



A( 7 ) = 



7-2- 1 

2r Ba 



■ sin 7 cos 7 



(24) 



(see their eqs. 1, 2, 4, and 16). In Purcell's (1979) heuristic derivation of the Barnett dissipation rate, he obtained the 
rate at which the rotational energy E decreases. In other words, he obtained the drift coefficient A(E). Since q oc E, 
there is no additional contribution to the drift coefficient arising from Ito's formula (|17[) when the Langevin equation 
is written in the variable q. However, the variable in Lazarian & Roberge (1997) is the angle 7, which is a non-linear 
function of E (eq. 3]). The additional contribution to the drift coefficient was not included in their analysis. Despite 
this error, the above general argument should still yield vanishing D and dD/d 7 at 7 = 7r/2, and thus no thermal 
flipping. 

In the vicinity of 7 = 7r/2 (corresponding to q = 1), the diffusion coefficient calculated by Lazarian & Roberge (1997; 
their eq. 18) is 



TBar-D(7) = ^^7 = ^/2) + 1 



(ra-l)fc- 



TBarC(7 = 7r/2) \ ( 7 - 7f/2) 



(25) 



If D(-f = 7r/2) is taken to be zero, then D oc (7 — 7r/2) 2 , as required. In this case, thermal flipping does not occur. 
Lazarian & Roberge (1997) argued that D should be smooth with respect to 7 and demanded that d 2 D jd^ 2 exist for 
all k. This condition is satisfied if D(-f = ir/2) = or if D oc k n with n < — 1. Lazarian & Roberge (1997) chose 
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D oc k 1 , which admits thermal flipping but is inconsistent with the requirement that D(j) falls off at least as quickly 
as (7 — 7r/2) 2 for 7 near ir/2. [There appear to be some typographical errors in Lazarian & Roberge 1997. In their 
eq. 19, D oc A: -1 / 2 rather than In their eq. 18, D{^ = tt/2) = 1 rather than falling off as fc -1 / 2 or fc -1 .] 

Lazarian & Roberge (1997) tested their result for the diffusion coefficient by numerically evolving their Langevin 
equation for a large number of Barnctt timescales and computing the average value of the internal alignment factor 



!x= 2 



(cos 2 7) - - 



(26) 



This can also be evaluated by simple integration for a thermal distribution (their eq. 10). The results of their 
simulations agreed to high accuracy with the direct calculations. They adopted the wrong Langevin equation but the 
correct thermal equilibrium distribution function. Their success with the test indicated that they solved equation (|10[) 
correctly given their drift coefficient, but this drift coefficient does not describe Barnett dissipation when angle 7 is 
taken as the variable. 

As a confidence-building check on the conclusion that thermal flipping is prohibited, I numerically evolved the 
Langevin equation for the case that k = 1 and r 2 — 1.5, taking C — 0. A fixed time step size is attempted at each 
step. Sometimes this results in overshooting q = 1; in these cases, smaller steps are tried until the resulting q exceeds 
1. These overshooting incidents become fractionally less common as the base step size is decreased (from 10 TBar 
to 10 _5 TBar)- The total duration of a simulation is about 10 5 TBar- At no time, for any of the base step sizes, did q 
ever overshoot r%. Incidentally, the simulations yielded the correct value for the alignment factor Qx, although the 
convergence was slower than for the simulations in Lazarian & Roberge (1997). 

The above argument that D and dD/dq both vanish at q = r 2 made no reference to the form of A(q). Thus, this 
conclusion also holds when more realistic Barnett dissipation rates are adopted, and even for grains lacking dynamic 
symmetry. In all of these cases, A(q = r%) =0, since the grain is in steady rotation when q = r 2 . Thus, q = r 2 is a 
natural boundary for the most general treatment of Barnett relaxation, if dD/dk exists for all k. Since a grain lacking 
dynamic symmetry must reach (and, in general, cross) q = r 2 in order to flip (see §2.5.2 of Weingartner & Draine 
2003), thermal flipping associated with Barnett relaxation is ruled out generally. The only caveat is that dD/dk must 
exist for all k. Although this seems natural, a detailed model of Barnett relaxation would be needed to confirm that 
this condition is indeed satisfied. 

The discussion here has focused on Barnett relaxation, since Barnett dissipation appears to dominate inelastic 
dissipation for most thermally rotating grains (Lazarian & Efroimsky 1999), especially when nuclear paramagnetism 
contributes. However, the argument against thermal flipping applies equally well for inelastic relaxation. Lazarian & 
Efroimsky (1999) did not constrain the form of the dissipation rate near q — r 2 , but Sharma et al. (2005) found the 
same form as in equation {2} for the special case of an oblate spheroid. 

3. CONCLUSION 

In conclusion, it appears that thermal flipping is not possible, so long as dD / dk exists for all k and the inertia tensor 
does not vary with time. A detailed model of Barnett relaxation is needed to examine the behavior of dD/dk. 

Because of grain vibrations, the inertia tensor exhibits continual, small variations. As a result, the location of the 
natural boundary at q = r% = wanders slightly (B.T. Draine, private communication). Further work is needed 

to examine whether this can give rise to flips and, if so, at what rate. 

External processes (e.g., gas atom impacts) may also induce flips (with accompanying changes in J), but have 
recently been neglected in comparision with internal relaxation (e.g., Weingartner & Draine 2003). These now merit 
further scrutiny as well. 
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